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We report a measurement of the time-dependent CP-asymmetry parameters 5 and C in color- 
suppressed B° — > D^*'°h° decays, where h° 'man , r\, or id meson, and the 75° decays to one of the 
CP eigenstates K + K~, K° s tt° , or ifgw. The data sample consists of 383 x 10 6 T(4S) -> SB decays 
collected with the BABAR detector at the PEP-II asymmetric-energy B factory at SLAC. The results 
are S = -0.56 ± 0.23 ± 0.05 and C = -0.23 ± 0.16 ± 0.04, where the first error is statistical and the 
second is systematic. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



Measurements of time-dependent CP asymmetries in 
B° meson decays, through the interference between de- 
cays with and without B°-B° mixing, have provided 
stringent tests on the mechanism of CP violation in the 
standard model (SM). The time-dependent CP asym- 
metry amplitude sin 2/3 has been measured with high 
precision in the b — ► ccs decay modes [![, where 
/3 = -a,rg(V c dV* b /VtdV t * b ) is a phase in the Cabibbo- 
Kobayashi-Maskawa (CKM) quark-mixing matrix Q ■ 

In this Letter, we present a measurement of the time- 
dependent CP asymmetry in B° meson decays to a neu- 
tral D meson and a light neutral meson through a b — » IJd 
color-suppressed tree amplitude. Interference between 
decay amplitudes with and without B°-B° mixing con- 
tribution occurs if the neutral D meson decays to a CP 
eigenstate. The measured time-dependent asymmetry is 
expected to be different from sin 2/3 measured in the char- 
monium modes due to the subleading amplitude b — > ucd, 
which has a different weak phase. This amplitude is sup- 
pressed by VubV* d /V c bV* d ~ 0.02 relative to the leading 
diagram. Therefore the deviation is expected to be small 
in the SM [llil. 

Many other decay modes that have significant contri- 
bution from loop diagrams have been studied [B| to con- 
strain or discover new physics due to unobserved heavy 
particles in the loop diagrams in B decays. This kind 
of new physics would not affect the decays presented in 
this Letter because only tree diagrams contribute to these 
modes. However, i?-parity- violating (fi p ) supersymmet- 
ric processes 0, 0] could enter at tree level in these de- 
cays, leading to a deviation from the SM prediction. 

The analysis uses a data sample of 348 fb _1 , which cor- 
responds to (383 ± 4) x 10 6 T(45) decays into BB pairs 
collected with the BABAR detector at the asymmetric- 
energy e + e~ PEP-II collider. The BABAR detector is de- 
scribed in detail elsewhere Q. We use the GEANT4 sim- 
ulation toolkit Q to simulate interactions of particles 
traversing the BABAR detector, and to take into account 
the varying detector conditions and beam backgrounds. 

We fully reconstruct B° mesons decaying into a 
CP eigenstate in the following channels: D^°tv° (D° 
K+K- 
D*° - 



, K^iu) [TO] and D^°r) {D° -> K + K~) with 
D 7r°, and D°uj (D° -> K+K~, K° s uj, K°ir°). 
From the remaining particles in the event, the vertex of 
the other B meson, B tag , is reconstructed and its flavor 
is identified (tagged). The proper decay time difference 
At = tcp — <tag, between the signal B (tcp) and B tag 



(t-t&g) is determined from the measured distance between 
the two B decay vertices projected onto the boost axis 
and the boost (/37 = 0.56) of the center-of-mass (cm.) 
system. The At distribution is given by: 



F±(At) 



-\At\/r 



4r 



[1 T Aw± 



(1) 



(1 - 2w)(r] f Ssm(AmAt) - C cos(AmAi))] , 

where the upper (lower) sign is for events with Stag be- 
ing identified as a B° (B°), iif is the CP eigenvalue of 
the final state, Am is the B°-B° mixing frequency, t is 
the mean lifetime of the neutral B meson, the mistag pa- 
rameter w is the probability of incorrectly identifying the 
flavor of Stag, and Aw is the difference of w for B° and 
B°. The neural- network based tagging algorithm [l2j has 
six mutually exclusive categories and a measured total 
effective tagging efficiency of (30.4 ± 0.3)%. Neglecting 
CKM-suppressed decay amplitudes, we expect the CP 
violating parameters S — — sin 2/3 and C = in the SM. 

The event selection criteria are determined by maxi- 
mizing the expected signal significance based on the sim- 
ulation of signal and generic decays of BB and e + e~ — > 
qq (q = u,d,s,c) continuum events. The selection re- 
quirements vary by mode due to different signal yields 
and background levels. 

A pair of energy clusters in the electromagnetic 
calorimeter (EMC), isolated from any charged tracks and 
with a lateral shower shape consistent with photons, is 
considered as a ir° candidate if both cluster energy de- 
posits exceed 30 MeV and the invariant mass of the pair 
is between 100 and 160MeV/c 2 . Charged tracks are con- 
sidered as pions, except for those used in D° — * K + K~ 
reconstruction, where the kaons must be consistent with 
the kaon hypothesis [HI]. We reconstruct 77 mesons in 
77 and 7T + 7r~7r modes. Each photon is required to have 
an energy exceeding 100 MeV and, when combined with 
any other photon in the event, to not have an invari- 
ant mass within 5MeV/c 2 of the n° nominal mass [l4j |. 
The invariant mass is required to be within approx- 
imately 30 MeV/ c 2 (8MeV/c 2 ) of the ?y nominal mass 
for 7] — > 77 (77 — > 7T + 7r~7r°). Both 7r° and 77 — > 77 
candidates are kinematically fitted with their invariant 
masses constrained at their respective nominal values. 
The u> — > 7T + 7r~7r° candidates are accepted if the in- 
variant mass is within approximately 22 MeV/ c 2 of the 
nominal u> mass, depending on the D° decay mode. The 
Kg — * 7r + 7r~ candidates are required to have an invari- 
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ant mass within 10MeV/c 2 of the K® nominal mass and 
X 2 probability of forming a common vertex greater than 
0.1%. The distance between the K® decay vertex and the 
primary interaction point projected on the plane perpen- 
dicular to the beam axis is required to be greater than 
twice its measurement uncertainty. 

The vector meson uj is fully polarized in D° — ► K®uj 
decays. Two angular distributions of the uj decay are 
used to discriminate against background: (a) cos 9®, de- 
fined in the uj rest frame, the cosine of the angle between 
the D° direction and the normal to the decay plane of 
uj — > 7r + 7T~7r , and (b) cos6>^, the cosine of the angle 
between the direction of one pion in the rest frame of 
the remaining pion pair and the direction of the pion 
pair. The signal are distributed according to cos 2 9® and 
1 — cos 2 Op, while the background distributions are nearly 
uniform. We require | cos6*^| > 0.4 and | cos 8p\ < 0.9. 

For the D° in D*° — > D°ir°, the invariant mass of the 
D° candidate is required to be within 30 MeV/c 2 of the 
world-average D° mass. For the D° in B° — > D°h°, the 
invariant mass window is tightened, ranging from ±14 
to ±29MeV/c 2 , depending on the mode. In both cases 
the D° is kinematically fitted with its mass constrained 
at its nominal value. The invariant mass difference be- 
tween D*° and D° candidates is required to be within 
±2.7MeV/c 2 of the nominal value. For B° -> D*°ir° 
with D° -> K° a uj, we require |cos#^| > 0.4, where 9* B 
is the angle between the momenta of the B° and the 7r° 
from the D*° in the D*° rest frame. 

The signal is characterized by the kinematic vari- 
ables m E s = \J (s/2 + p • Vb) 2 / e I ~ Ps> and AP = 
Eg — E^ eaxn , where the asterisk denotes the values eval- 
uated in the cm. frame, the subscripts 0, beam and 
B denote the e + e~ system, the beam and the B can- 
didate, respectively, and y/s is the cm. energy. We re- 
quire toes > 5.23 GeV/c 2 . The AE distribution for signal 
events is asymmetric and varies by decay mode. Depend- 
ing on the mode, the lower (upper) boundary of the AE 
selection window varies from —95 to — 35MeV (+35 to 
+85MeV). The reconstructed |Ai| and its uncertainty 
CAt are required to satisfy |Ai| < 15 ps and cta* < 2.5 ps. 

The background from continuum qq production is sup- 
pressed based on the event topology. In the cm. frame, 
the B mesons are produced nearly at rest and decay 
isotropically, while the quarks in the process e + e~ — > qq 
are produced with large relative momentum and result in 
a jetlike topology. The ratio of the second to zeroth order 
Fox- Wolfram moments [IH , determined from all charged 
tracks and clusters in the EMC with energy greater than 
30MeV, must be less than 0.5. The qq background 
is further suppressed by a Fisher discriminant T [lj|, 
constructed with the following variables, evaluated in 
the cm. frame: (a) L2/L0 where Li — YljPjl c os9* | 4 , 
summed over the remaining particles in the event after 
removing the daughter particles from the B°, p* is the 
momentum of particle j and 9* is the angle of the momen- 



tum with respect to the B° thrust axis [l7| ; (b) | cos 9^ | , 
where 9^ is the angle between the B° thrust axis and the 
thrust axis of the rest of the event; (c) | cos 2 9* B \, where 
9* B is the angle between the beam direction and the di- 
rection of the J5°; (d) total event thrust magnitude; and 
(e) total event sphericity [l8l ]. 

For B° — > D°oj decays, we add two angular variables 
to T\ cos 9^ and cos 9^ , analogous to cos 9^ and cos 9^ 
in D° — > KgOJ. The signal distributions for the B° sys- 
tem are the same as those in the D° system. The back- 
ground distributions are close to 2 — cos 2 9^ and uniform 
in cos6>£. The requirement on T depends on the back- 
ground level in each mode; the signal selection (back- 
ground rejection) efficiency is 60%-86% (72%-94%). 

Within each reconstructed decay chain, the fraction of 
events that have more than one candidate ranges from 
less than 1% to about 10%, depending on the mode. We 
select one candidate with the most signal-like Fisher dis- 
criminant value for each mode. A total of 1128 events 
are selected, of which 751 are tagged (the absolute value 
of the flavor-tagging neural-network output greater than 
10% of the maximum). 

The signal and background yields are determined by 
a fit to the uies distribution using a Gaussian distribu- 
tion for the signal peak and a threshold function [3] 
for the combinatorial background. We obtain 340 ± 32 
signal events (259 ± 27 tagged). The contribution from 
each mode is shown in Table HI and the toes distribu- 
tions are shown in Fig. [TJ We investigate potential back- 
grounds that might peak in the toes signal region by 
studying data in the D° mass sideband (outside a window 
of ±3 standard deviations of the mass peak) and simu- 
lated e+e~ -> BB events. We estimate that (0.8 ± 2.6)% 
of the CP-even signal yield and (5.4 ± 2.2)% of the CP- 
odd signal yield are background, based on the simulation. 
Approximately half of the peaking background found in 
simulation is from B~ — > D°p~(— > 7r°7r~) with a low 
momentum tt~ . Other sources include B° — > ir + ir~ir° 
and B° — > D^*>°h°, with D° decaying to a flavor eigen- 
state, e.g., K~ir + . We find that the peaking background 
from the D° mass sideband data in CP-even modes is 
consistent with the simulation. For CP-odd modes, we 
find a larger peaking component in D° sideband data 
than expected from simulation. Therefore we increase 
the estimated total peaking background fraction for CP- 
odd events to (11 ± 6)% to account for the excess found 
in the D° sideband data. We estimate that 65% of the 
peaking background arises from charmless decays with 
potentially large CP-violating asymmetries. We account 
for this possibility in the systematic uncertainty. 

In order to extract CP violating parameters S and 
C, we fit the toes and At distributions of the 751 
tagged events using a two-dimensional probability den- 
sity function (PDF) that contains three components: sig- 
nal, peaking background and combinatorial background. 
The Toes distribution is described in the previous para- 
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TABLE I: Signal yields. Uncertainties are statistical only. 
The CP parity of the D is indicated in the column of Dqp- 
The combined value is from a simultaneous fit to all modes. 



U 


+ 1 (CP 


even) 




-1 (CP odd) 


Mode 


D C p 


-^signal 


Mode 


D C p 


N signal 






26.2 ±6.3 


Dkrk 


± 


104 ± 17 






40.0 ±8.0 


DkkVjj 


+ 


28.9 ±6.5 






23.2 ±6.8 


D°KKV3n 


+ 


14.2 ±4.7 


+ 


23.2 ±6.3 


D° kk lo 


+ 


51.2 ± 8.5 




+ 


9.8 ±3.5 


D K°w* 




5.5 ±3.3 




+ 


6.8 ±2.9 








Combined 




131 ± 16 






209 ± 23 


Total 










340 ± 32 



graph. Its parameters are free in the fit. The peaking 
background is assumed to have the same toes shape as 
the signal. The signal decay-rate distribution shown in 
Eq. |T]) accounts for dilution due to an incorrect assign- 
ment of the flavor of -Btag, and is convolved with a sum 
of three Gaussian distributions, parameterizing the core, 
tail and outlier parts of the At resolution function [l3j ]. 
The widths and biases of the core and tail Gaussians are 
scaled by a At- The biases are nonzero to account for 
the charm meson flight from the -Btag vertex. The out- 
lier Gaussian has a fixed mean (0 ps) and width (8 ps) 
to account for poorly-reconstructed decay vertices. The 
mistag parameters and the resolution function are deter- 
mined from a large data control sample of B° — > D^*'h + 
decays, where h + is a 7r + , p + , or meson. The B° life- 
time and mixing frequency are taken from Q . 

We use an exponential decay to model the At PDF 
of the peaking background. We account for possible CP 
asymmetries in the systematic uncertainty. The At PDF 
for combinatorial background consists of a term with zero 
lifetime to account for the qq contribution, and an oscil- 
latory term whose effective lifetime and oscillatory co- 
efficients are free parameters in the fit to account for 
possible CP asymmetry in the background. The sum of 
a core Gaussian and an outlier Gaussian is sufficient to 
model the resolution function. The combinatorial back- 
ground parameters are determined predominately by the 
events in the Toes sideband. The final PDF has 25 free 
parameters for fitting to all modes and tagging categories 
simultaneously. 

We obtain S = -0.56 ± 0.23 ± 0.05 and C = -0.23 ± 
0.16 ± 0.04, where the first errors are statistical and the 
second are systematic. The statistical correlation be- 
tween S and C is p = —2.4%. The At distribution 
projections and the asymmetry (A — [N B o tag (At) — 

N S ot ag (^)]/[NBH, g (At) ± AWg( A *)D for the events 
in the signal region are shown in Fig. O We check 
the consistency between CP-even and CP-odd modes 
by fitting them separately and find (statistical errors 
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FIG. 1: The toes distributions with a fit to (a) the CP-even 
and (b) the CP-odd modes combined in the data. The solid 
curve represents the overall PDF projection and the dashed 
curve represents the background. 




5 10 
At (ps) 



5 10 
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FIG. 2: The At distributions and asymmetries for (a,b) CP- 
even and (c,d) CP-odd events in the signal region (toes > 
5.27GeV/c 2 ). In (a) and (c), the solid points with error 
bars and solid curve (open circles with error bars and dashed 
curve) are P°-tagged (_B°-tagged) data points and At projec- 
tion curves. Shaded areas (P°-tagged) and the dotted lines 
(P°-tagged) are background distributions. In (b) and (d), the 
solid curve represents the combined fit result, and the dashed 
curve represents the result of the fits to CP-even and CP-odd 
modes separately. 



only) Seven = -0.17 ± 0.37, 5 odd = -0.82 ± 0.28, and 
Ce V en = -0.21 ± 0.25, C od d = -0.21±0.21. The difference 
between iS OV en an d <? dd ^ s 0-65 ± 0.46, less than 1.5 stan- 
dard deviation from the expected value, zero. We also 
find that the differences between h° — * 77 and h° — * tttttt 
modes are less than 0.1 in C and S. 

The SM corrections due to the sub-leading-order dia- 
grams are different for Dcp+ and Dcp- [3]. Therefore, 
we also perform a fit allowing different CP asymmetries 
for L>cp + and D C p— We obtain S+ = -0.65±0.26±0.06, 
C+ = -0.33 ± 0.19 ± 0.04, p+ = 4.5%, and S_ = 
-0.46±0.45±0.13,C_ = -0.03±0.28±0.07, p_ = -14%. 

The dominant systematic uncertainties are from the 
peaking background and the toes peak shape uncertain- 
ties (0.04 in S and 0.03 in C). For the former, we vary 
the amount of the peaking background according to its 
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estimated uncertainty, and vary the CP asymmetry of 
the charmless component between ± sin 2/3 of the world- 
average value. We study the latter effect using an alter- 
native line shape [20| taking into account a possible non- 
Gaussian tail in the toes distribution. Other systematic 
uncertainties typically do not exceed 0.01 in S or C, and 
come from the following sources: the assumed parame- 
terization of the At resolution function; the uncertainties 
of the peaking background; toes width and the combi- 
natorial background threshold function; B° lifetime and 
mixing frequency; the beam-spot position; and the in- 
terference between the CKM-suppressed b — * Tied and 
CKM-favored b — > cud amplitudes in some B teig final 
states, which gives deviations from the standard time 
evolution function Eq. ([I]) [2l[. Uncertainties due to the 
vertex tracker length scale and alignment are negligible. 
Summing over all systematic uncertainties in quadrature, 
we obtain 0.05 for S and 0.04 for C. 

In conclusion, we have measured the time-dependent 
CP asymmetry parameters S = —0.56 ± 0.23 ± 0.05 and 
C = -0.23 ± 0.16 ± 0.04 from a sample of 340 ± 32 B° -> 
Dgph signal events. The result is 2.3 standard devia- 
tions from the CP-conserving hypothesis S = C = 0. The 
parameters S and C are consistent with the SM expecta- 
tion, i.e., the world average — sin 2/3 = —0.725 ±0.037 [f| 
and zero, respectively. 
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